Superconductivity in heavily vacant diamond 
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Using first principle electronic structure calculations we investigated the role of substitutional 
doping of B,N,P,A1 and vacancies (V) in diamond (XaCi_a). In the heavy doping regime, at about 
~ 1 — 6% doping an impurity band appears in the mid gap. Increasing further the concentration 
of the impurity substitution fills in the gap of the diamond host. Our first principle calculation 
indicates that in the case of vacancies, a clear single-band picture can be employed to write down an 
effective one band microscopic Hamiltonian, which can be used to further study various many-body 
and disorder effects in impurity band (super)conductors. 
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I. INTRODUCTION 

Diamond has a number of unique attributes that make 
it highly suited as a gem stone. It is the hardest known 
material which can only be scratched by another dia- 
mond. The thermal conductivity of diamond is the high- 
est among all materials [l|, [2| • Irradiation of diamond by 
various particles (e.g. electrons, neutrons, a particles) 
followed by annealing to repair damaged sj? bonds gives 
rise to fascinating colours of diamond, which are due to 
the so called colour centers [3]]. 

Diamond is also a material with semiconductor prop- 
erties that are superior to Si, Ge, or GaAs, which are now 
commonly used. The use of diamond in electronic appli- 
cations is not a new idea, but limitations in size and con- 
trol of properties restricted the use of diamond to a few 
specialized applications. The vapor-phase synthesis of 
diamond, however, has facilitated serious interest in the 
development of diamond-based electronic devices. The 
process allows diamond films to be laid down over large 
areas. Both intrinsic and doped diamond films have a 
unique combination of extreme properties for high speed, 
high power and high temperature applications [2, |j, l5| . 

Ekimov and coworkers _7] and subsequently Takano ei. 
al. Q used chemical vapor deposition (CVD) to synthe- 
sized B-doped diamond. Doping diamond by low concen- 
tration of typically 10^^ — 10^^ cm~'^ boron atoms gives 
rise to acceptor level, rendering it to a p type semiconduc- 
tor 0]. Increasing the doping level to [B]/[C]> 5000 ppm 
in the gas phase induces metallic conductivity in dia- 
mond [^ . Further increasing the doping rate to the scale 
of n > 10^^ cm~^, i.e. ~ few %, makes it supercon- 
duct at low temperatures [3, Q . Increasing the doping 
rate amounts to bringing the boron atoms closer, and al- 
low them to overlap more effectively, which broadens the 
acceptor levels in to decent bands of electrons [Ifl, El , 
which are resp onsible for metallic and superconducting 
properties [11 [H, [H [H, [H [H, [I^ . 

Also some authors have investigated the doping of sili- 
con with boron, aluminum and phosphorus [18lll9|. Ex- 
periment has showed that the transition temperature for 
B-dopcd Si is T^ ~ 0.35 K. Therefore doping C with 
boron gives more than an order of magnitude larger 
Tc ~ 4 K compared to B doped Si. Hence we choose 



the host material to be carbon rather than any other 
element in the same group. 

Now the question arises, what other elements can be 
doped into diamond in the regime of heavy (~ 1 — 10%) 
doping which can possibly lead to higher superconduct- 
ing Tc- In the example of high temperature cuprate su- 
perconductors (HTSC) [20] , an effective one-band model 
for the so called Zhang-Rice singlet [21j| can be written 
down in terms of the hole states of the 02p and CuSd. 
In Xa_C\-a case also the effective impurity band has a 
mixed Xnp and C2p character. Here n = 2 for X=B,N 
and n — Z for X=A1,P. In the case of X=V the pic- 
ture is even simpler. The metallic band in the middle 
of the diamond gap well isolated from both bands is al- 
most entirely due to the nearest neighbor C atoms sur- 
rounding the vacant site. In this work we present our 
preliminary result on comparison of the impurity band 
formation for the above elements for various doping rates 
a = 1/128,1/54,1/16. 



II. METHOD OF CALCULATION 



In this study we used the Plane Wave-Pseudopotential 
Quantum-ESPRESSO code ^. We used Density Func- 
tional Theory with General Gradient Approximation 
(GGA). The GGA exchange-correlation functional which 
has been used is PBE [2>|]. We employed ultrasoft pseu- 
dopotential [2J] to describe electron-ion interaction. The 
energy cutoff for expansion of wave function in plane 
wave was 25 Ry, and 150 Ry has been used for expan- 
sion of charge density. We chose 2x2x2,3x3x3 
and 4x4x4 supercell with one vacancy or with a de- 
fect, that is XttCi_„ with a = 1/128,1/54,1/16. B, N, 
Al, P are chosen to substitutionally replace one carbon 
atom in a diamond structure. The k-point is sampled 
according to table ID To accelerate electronic structure 
calculation, we use Methfessel and Paxtons Fermi-level 
smearing method (width 0.01 Ry) [25]. The impurity has 
been only substituted by one of carbon atom in diamond 
structure without any relaxation. 



TABLE I: The k-point sampling for different supercells. 



Supercell 


k-point grid 


number of 
atomic sites 


2x2x2 


6x6x6 


16 


3x3x3 


4x4x4 


54 


4x4x4 


3x3x3 


128 




FIG. 1: The density of state (DOS) for various concentrations 
of B, AI, N, P and V impurities in diamond. 



III. RESULTS AND DISCUSSION 

The results of this work has been summarized in Fig. [T] 
where we have plotted the density of state (DOS) for 
various impurities and different impurity concentrations. 
We have shifted the data such that the Fermi level in 
all cases is at Ep — 0. The impurity band formation for 
the vacancy and Aluminum are the most clear among the 
cases studied here. Their bandwidths are between 1.5 eV 
to 0.5 eV for different impurity concentration. 

For the perfect diamond lattice the DOS for all sizes 
of the unit cell coincide and give the clean diamond gap. 
However substitutional doping with B and Al starts to 
create acceptor levels on top of the valence band which 
broaden into bands in the regime of few percent impurity 
concentration. Also in this regime the relevant orbitals of 



the impurity atoms will hybridize with the C2p orbitals, 
and the impurity band has a mixed character, similar 
to the case of Zhang-Rice singlet of cuprates [2l|. It is 
clearly seen that the impurity band peak in the case of 
Al is stronger than B. Also the value of the DOS at the 
Fermi level (po) is larger for Al doping than B doping. N 
and P similarly create donor levels at the buttom of the 
conduction band. Again qualitatively one can see that P 
tends to give a sharper impurity band peak than the N. 
Here also the impurity band arises from a combination 
of nearest neighbor C2p and N2p or P3p bands. 

In the case of vacancy, the story is different. First of all, 
since V is neither acceptor, nor donor level, the resulting 
band will be in the middle of the gap, well isolated from 
the valence and conduction band. Secondly the impu- 
rity band arises from the n.n. carbon atoms surrounding 
the vacancy. This qualitative picture can be inferred by 
looking at the orbital and site resolved partial DOS (not 
shown here). Therefore for the effective one band model 



of the impurities proposed by Baskaran [ll[ , doping by 
V seems to be more suitable than the other cases studied 
in this investigation. 

Note that in this study we have ignored issues like 
the formation of various complexes. For example, ni- 
trogen doping in diamond usually leads to the formation 
of nitrogen- vacancy complex [26|. Ignoring such com- 
plications, one can use the argument of Mott to get a 
rough estimate of the typical density needed to make the 
the resulting half-filled impurity band superconduct: The 
critical concentration needed for metalization is given by 
a^ric = 1/4, where the Bohr radius ub of the impurity 
can be estimated from the binding energy Eb of the im- 

2 

purity levels as Ub = 2^^ e ■ -'^°'' ^ypi^ally 0.5 — 1.5 eV 
binding energies of say, Al, V, the Bohr radius will be 1—3 
atomic units, giving a critical concentration on the scale 
of ~ 10^^ cm~^ or 5 — 15 percents. For the concentrations 
affordable in our calculations, at a = 1/16 « 6% non of 
the impurities studied here fills in the gap. However Cu 
doping at nearly 6% already metalizes the diamond [23] . 
According to a disordered RVB mechanism suggested 
by Baskaran [ll| , for such an effective single band at half 
filling the critical temperature for the superconductivity 

is given by UbTc « -^epo-' , where J is superexchange, 
W the band width, and po is the DOS at the Fermi level 
obtained from the DFT bands. The effective one band 
model must have large enough W to escape the Mott- 
Hubbard splitting in the large U limit, such that the 
half-filled band picture remains valid. Also the disor- 
dered nature of the impurity centers will start to localize 
the states close to the impurity band edges. Again the 
bandwidth W must be wide enough such that the mobil- 
ity edge will not cross the Fermi level. 

If one assumes that the superconducting mechanism 
remains the same for various impurities X studied here, 
vacancy and Al doping offer a more clear picture of the 
half-filled band undergoing Anderson-Mott to RVB su- 
perconducting scenario of Baskaran [11|. The enhance- 
ment of po observed for the case of similar concentration 



of Al and V is advantageous in giving larger transition 
temperature than the case of Boron. Note that larger po 
even within the BCS picture is an advantage of V and Al 
doping compared to doping by B. 

In terms of practical fabrication, utilizing the heavily 
vacant diamond may offer a new method in addition to 
high temperature high pressure techniques used in pro- 
duction of CVD diamond doped with various elements. 
Irradiation by different particles may offer a method of 
producing high concentration of vacancies in diamond 
which can possibly lead to higher transition temperature 



than the CVD diamond doped with elements. 
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